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O ! Abstract 

cn 
o 

^ I The angular power spectrum of the Sunyaev-Zeldovich (SZ) effect is calculated 

Q^' in the ACDM cosmological model with the aim of investigating its detailed depen- 

Q . dence on the cluster population, gas morphology, and gas evolution. We calculate 

|H I the power spectrum for three different mass functions, compute it within the frame- 

^ ■ work of isothermal and polytropic gas distributions, and explore the effect of gas 

evolution on the magnitude and shape of the power spectrum. We show that it is 
indeed possible to explain the 'excess' power measured by the CBI experiment on 
^ ' small angular scales as originating from the SZ effect without (arbitrary) rescaling 

^ . the value of erg, the mass variance parameter. The need for a self-consistent choice 

of the basic parameters characterizing the cluster population is emphasized. In par- 
ticular, we stress the need for a consistent choice of the value of ug extracted from 
fitting theoretical models for the mass function to the observed cluster X-ray tem- 
perature function, such that it agrees with the mass-temperature relation used to 
evaluate the cluster Comptonization parameter. Our treatment includes the explicit 
spectral dependence of the thermal component of the effect, which we calculate at 
various frequencies. We find appreciable differences between the nonrelativistic and 
relativistic predictions for the power spectrum even for this superposed contribution 
from clusters at the full range of gas temperatures. 
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1 Introduction 



The Sunyaev-Zeldovich (SZ) effect - the change in the cosmic microwave back- 
ground (CMB) intensity that results from Compton scattering in a cluster of 
galaxies - is a unique cosmological probe which is essentially independent of 
the cluster redshift. The effect has already been imaged in some ~ 60 clusters, 
yielding important information on such quantities as the gas mass fraction 
in clusters, and the Hubble constant, Hq. Multi-frequency measurements of 
the effect in many more clusters are expected in the near future when new 
SZ projects - both ground-based and stratospheric - will be operational. In 
addition to the detection of the effect in a large number of clusters, planned 
sky surveys, particularly by the Planck satellite, will likely map the CMB 
anisotropy induced by the effect. The much larger SZ database and the mea- 
surement of this anisotropy will greatly expand the scope of the effect as a 
cosmological probe due to its strong dependence on the cosmological and large 
scale parameters, and the morphology and evolution of clusters. 

Extensive work in the last two decades (since it was first modeled [Rephaeli 
1981] in the context of a simple model for IC gas evolution) has established 
the main features of SZ anisotropy in various cosmological, dark matter (DM), 
and cluster evolution models. For example, the SZ power spectrum was calcu- 
lated in fiat cold dark matter {e.g. , SCDM, Atrio-Barandela & Miicket 1999), 
ACDM {e.g., Komatsu & Kitayama 1999), and open low-density (OCDM, 
Molnar Sz Birkinshaw 2000) models (for reviews, see Rephaeh 1995a, Birkin- 
shaw 1999, and Carlstrom et al. 2002). However, results for the power spectrum 
of this anisotropy and cluster number counts differ considerably even when cal- 
culated for the same cosmological and DM models, mostly due to differences 
in the modeling of the cluster mass function and gas properties, but in some 
cases also due to inconsistent choice of some of the many free parameters. 
There is a need to clarify some of these differences as part of the basic goal 
of a quantitative and realistic treatment of the SZ-induced anisotropy. This is 
clearly desirable also in order to optimize strategies for the spectral and spa- 
tial mapping of the anisotropy in the many surveys planned by ground-based 
and satellite experiments. 

The measured values of CMB intensity differences across the sky naturally 
consist of contributions from all sources of anisotropy - primary, secondary, 
and tertiary. Most investigated is, of course, the primary anisotropy, orig- 
inating in the very early universe, for which the SZ contribution could be 
regarded as a foreground 'contamination'. The latter, however, is very im- 
portant for better understanding of the post-recombination universe, and - 
in particular - to the global distribution of galaxy clusters and their physi- 
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cal properties. The primary anisotropy has to be precisely known before the 
SZ (and other contributions) can be determined on angular scales of a few 
arcminutes. This is particularly important now that an 'excess' power was 
measured by the CBI experiment (Mason et al. 2002), at an estimated 3. la- 
level, over what is predicted in standard models for the primary anisotropy at 
multipoles ^ = 2000-3500. This 'excess' can possibly be due to the integrated 
SZ effect (Bond et al. 2002). The power spectrum was recently determined 
also by measurements with ACBAR (Goldstein et al. 2002) at ~ 150 GHz. 
Measurements at this and higher frequencies - well beyond the CBI 26-36 
GHz band - are essential in order to exploit the unique spectral signature of 
the SZ effect as a powerful diagnostic tool by which it can be separated out 
from the primary anisotropy {e.g. , Rephaeli 1995a, Cooray et al. 2000). 

The calculation of the SZ power spectrum necessitates detailed modeling of 
the cluster mass function, intracluster (IC) gas properties and evolution. Con- 
sequently, the values of many parameters have to be specified in order to fully 
characterize clusters and their gaseous contents. Some of these parameters may 
be obtained from theoretical considerations (e.g. the critical overdensity for 
spherical collapse), while others are extracted from observations (e.^. , cluster 
X-ray temperatures), or by comparing theoretical models with observational 
data, yielding important normalization factors, such as a^. This important 
parameter is the rms density field smoothed over a scale 8/i~^ Mpc, where h 
is the Hubble constant in units of 100 km/(s Mpc). In particular, the cluster 
mass-temperature relation and cxg, which is usually estimated by means of the 
cluster X-ray temperature function, are coupled quantities which have to be 
consistently selected. 

Since the angular power spectrum of the SZ effect has already been calculated 
in numerous papers, we focus here on various aspects of the calculation which 
have not yet been sufficiently explored. These include explicit form for the 
temporal of IC gas, and a more general polytropic spatial distribution of its 
temperature. The description of the cluster population is also extended to in- 
clude two mass functions incorporating non-spherical collapse, in addition to 
the Press & Schechter mass function. We also discuss the need for a consistent 
choice of normalization of the cluster density and mass-temperature relation. 
In our claculation we include the full spectral dependence of the thermal effect 
which we calculate at five frequencies (using the exact relativistic treatment) 
below, near, and above the crossover frequency (where the thermal component 
vanishes). These calculations underline the importance of using the exact rel- 
ativistic expression for the change of intensity due to the SZ effect (Rephaeh 
1995b), especially so near the crossover frequency. 

Our basic computational approach is essentially similar to that presented in a 
few other papers {e.g. , Colafrancesco et al. 1997, Molnar & Birkinshaw 2000); 
thus, our discussion of most aspects of the methodology will be brief. Rather, 



3 



we will concentrate on issues that have attracted little or no attention. The 
paper is arranged as follows: in § 2 we lay out the computational method, 
specifying the different parameters and scalings required in the calculation of 
SZ anisotropy. Results for the angular power spectrum are presented in § 3, 
and discussed in § 4. 



2 Methodology 

The calculation of CMB anisotropy induced by the SZ effect requires modeling 
of IC gas, the cluster population, and their evolution, as well as the parameters 
of the background cosmological model. We begin with a brief description of the 
effect itself: The thermal component of the effect is described by the change of 
the spectral intensity of the CMB across the cluster, A/j. In the nonrelativis- 
tic approximation, Alt = ioUgoix), where = 2{kT)^/{hcf, T is the CMB 
temperature, y = J{kTe/'mc'^)naTdl is the Comptonization parameter, and 
go{x) is the well known (Sunyaev & Zeldovich 1972) spectral function. Here n 
and Te are the electron number density and temperature, ax is the Thomson 
scattering cross section, d£ is a line element along a line of sight (los) to the 
cluster, and x = hu/kT is the non-dimensional frequency. This nonrelativis- 
tic expression for A/^ is not sufficiently accurate at high frequencies and at 
gas temperatures higher than a few keV; a more exact relativistic description 
is then required (Rephaeli 1995b). As shown in the latter paper, the more 
exact relativistic calculation A/^ differs appreciably from the non-relativistic 
expression, particularly so near the crossover frequency, Xq. This has practical 
ramifications since xq is the optimal frequency for measurement of the cluster 
velocity (in the CMB frame) along the los, and the fact that this frequency is 
extremely suitable for CMB observations, as it allows the subtraction of (the 
dominant) thermal component from the full CMB signal. 

2.1 IC Gas Properties 

In all calculations it is assumed that IC gas is in a state of hydrostatic equilib- 
rium in the gravitational potential of the cluster. The density profile adopted 
here is the commonly used isothermal (3 model 



where Vf. is the core radius. For simplicity we take the cluster radial extent - 
roughly, the virialization radius - to be pr^ = lOrc. The core radius is scaled 
such that Tc — 0.15 Mpc for a local cluster with mass 10^^ Mq. 
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The gas mass is assumed to constitute a fraction / of the total cluster mass. 
For local clusters of mass 10^'^ MQ^h^^. wc will scale to the nominal value 
/ = 0.1, in accord with recent SZ results (Carlstrom et al. 2002). The fact 
that the gas is appreciably metal-enriched, with typical iron abundance of 
~ 0.2 — 0.3 solar, is indicative of galactic and intergalactic origins and varying 
relative contributions. Temporal evolution of the gas must then be taken into 
account; we do so by the parameterization f{M, z) oc M'H^{z), with t denot- 
ing the age of the universe at redshift z. Some evidence for the form of this 
evolution was seen in cluster measurements made during the Einstein Medium 
Sensitivity Survey (EMSS; Gioia et al. 1990, Henry et al. 1992). Analysis of 
the EMSS results led Colafranccsco & Vittorio (1994) to the deduction of the 
approximate values ^ = 1.45, and r] = 0.2, which we adopt here. Clearly, an 
increase of IC gas fraction with mass, and a decrease with z is indicated by 
these values. 

We consider general polytropic models Tg oc rf~^ which include the special 
case of isothermal distribution (7 = 1). The integrated y parameter along a 
los forming an angle 6 with a los pointing at the cluster center is given by 
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Here yo = ?^ii„(M,2)r„(M,z)r„(M,2), n = slf - <S = 3/37/2, 
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and 2-^1 is the Hyper geometric function. In the simpler isothermal case we 
adopt /3 = 2/3, and obtain 
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2.2 Cluster SZ Profile 



The calculation of the SZ angular power spectrum involves a number of stages. 
Here we only outline the method; further details can be found in previous 
works {e.g. , Colafrancesco et al. 1997, Molnar & Birkinshaw 2000). The rel- 
ative temperature change along a los to the cluster is 
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where x — hu/kT is the non-dimensional frequency, and y{6) is the Comp- 
tonization parameter along a los displaced an angle 6 from the center of the 
cluster. For isothermal gas with a King profile we have 



y{0) - 2-^ niO){z, M) T{0){z, M) r,{z, M) 



tan 



1 - {0/pOc) 



(5) 



where Oc denotes the angular diameter subtended by the cluster core. Ex- 
pressions for the temperature, density and cluster core radius are given in 
Colafrancesco et al. (1994) for flat cosmological models, equations (6), (9) and 
(10), and the equivalent relations in equations (6) and (7) in Colafrancesco 
et al. (1997) for low-density models. For polytropic gas the corresponding 
expression is 
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The spectral dependence of (the thermal component of) AT, which is given 
by the function s{x) — xcoth(a;/2) — 4, is only valid at low cluster tem- 
peratures. As noted, a relativistically correct expression has to be used for 
most of the cluster temperature range and at high frequencies. Approximate 
analytic expressions for the relativistic case have been obtained, applicable 
at temperatures ksT < IbkeV {e.g., Itoh et al. 1998, Shimon & Rephaeli 
2003). The computation of the SZ angular power spectrum involves an in- 
tegration over the mass-redshift distribution of the cluster population, and 
thus clusters in a wide range of temperatures - determined by the mass- 
temperature (and redshift) scaling - contribute to the power spectrum. The 
predicted range of gas temperatures can be estimated by using the usual scal- 



ing r= 7. 76/3" 
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[l + z) keV, in a flat = 1 CDM model 



(Molnar & Birkinshaw 2000). Since j3 assumes values in the range ~ 0.5 — 1, 
and M can be as high as 10^^ h^^ Mq, it follows that this relation yields high 
temperatures (> 30 keV) at z > 1. While a naive use of this scaling relation 
at ^ > 1 is likely to be improper - due (among other considerations) to fast 
radiative cooling in the high density cores of very massive clusters, if these 
existed at all already at such early times (contrary to expectations of hierar- 
chical growth by merging of sub-clusters) - it is clear that we have to use the 
relativistically correct spectral distribution of the SZ effect. It is nonetheless 
obvious that the mass function yields relatively low number of clusters with 
high temperatures, and that their overall contribution to the power spectrum 
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is likely to be small. 



2.3 Cluster Mass Functions 



We adopt the PS mass function (Press & Schechter 1974) as the 'standard' 
model for characterizing the cluster population. However reasonable and suc- 
cessful this function is, it still suffers from a number of drawbacks. From a 
theoretical point of view, the function has been derived using the oversimpli- 
fied assumption of spherical infall, which lies on weak theoretical grounds. Its 
failure to deal with secondary infall is yet another well-known shortcoming. 
It has been shown that other theoretically derived mass functions are more 
consistent with results of N-body simulations. Specifically, Lee & Shandarin 
(1999), and Sheth & Tormen (1999) (hereafter LS & ST, respectively) have 
derived mass functions incorporating non-spherical collapse. We calculate the 
SZ angular power spectrum based on all these three mass functions and assess 
the differences between the respective values of Q. 

For the mass function to correctly quantify the cluster abundance in the mass- 
redshift space, the normalization of the power spectrum of the density fiuc- 
tuation field must be determined. This is usually specified in terms of erg, 
the rms density field smoothed over a scale 8h~^ Mpc, a free parameter in 
the analytic mass function expression whose best- fit value can be extracted 
from the observed X-ray temperature function. The temperature function de- 
scribes the abundance of clusters N{T) with temperature T within a small 
interval dT; this is the differential form of the function. The integral form 
A^(> T) specifies the population of clusters with temperatures equal to or 
higher than T. Fitting a mass function to the temperature function requires a 
mass-temperature calibration, ascribing a temperature T{M) to a cluster with 
mass M. The parameter as obviously depends on the adopted calibration. We 
first explain this dependence intuitively, then point at likely consequences of 
using a normalization which stems from an inconsistent calibration. 

Consider two temperature-mass relations (TMR) denoted by Ti(M) and T2(M), 
satisfying Ti(M) > T2(M) for all M. Since any TMR predicts higher tempera- 
tures with increasing mass, i.e., the temperature is a monotonically increasing 
function of mass, the inverse mass-temperature relations (MTR), Mi(T) and 
M2(r), must satisfy Mi{T) < M2(T) for all T. Consider now an observed 
temperature function N{T)dT, which has to be related to the mass function 
(MF). The latter is characterized by the normalization factor cg, and a fit of 
the presumed MF to the observed temperature function using a TMR may be 
performed with the normalization erg as the fit parameter. The temperature 
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function may be written in terms of the MF as follows: 

N{T)dT^N[M{T)]—dM, (7) 



and 



N[M{T)] dM = N{T) ^dT^ N{T) dT ^, (8) 



where the last equality follows from the relation T oc M^/^. Given two MTRs 
the corresponding MF predictions are 



N[Mi{T)] dM = N{T) dT ^^^^^^ 



N[M2{T)] dM = N{T) dT 



2T 

(9) 

3M2(r) 

2T 



However, to be able to compare their predictions, the arguments of these 
functions must be identical. Defining Mi(T) = M2(T*) we obtain the following 
set of equations: 



n[Mi(r)] dM = n{T) dT ^^^^^^ 



n[M2(T*)] dM = n{T*) dT 



2T 

(10) 

3M2(r*) 

2T* 



Now the arguments of the two MFs are identical but the MFs predictions 
themselves are different, and since Mi(T) < Af2 (T)VT, it turns out that 
in order for the condition Mi(T) = M2{T*) to be satisfied, one must have 
T* < T. This means that the lower expression in equation set (10) predicts 
higher cluster abundances. We may thus conclude that a MF derived using a 
MTR which associates a lower mass with a given temperature would predict 
a higher cluster abundance at the corresponding mass. 

The MF is to a good approximation a monotonically increasing function of 
(Tg, technically due to the fact that ag appears in the denominator of the ex- 
ponential part of N{M). Only in the small mass range does the multiplicative 
term 1/(Ts dominates its behavior of the MF; in this range the MF is a de- 
creasing function of as [N{M) oc (Tg^exp{—6'^/2(T'^ajDl)]. Specifically, N{M) 
is an increasing function of as for ctm ■ < 5c, where ctm, and 5c are the 
(unnormalized) mass variance smoothed over mass scale M, the linear growth 
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factor of density fluctuations, and the critical density for spherical collapse, 
respectively. We can now predict the consequences of employing a normaliza- 
tion factor as extracted using a specific MTR, but using a different MTR for 
specifying the temperature ascribed to a cluster of mass M, required for the 
evaluation of the cluster y-parameter oc Tq. Note that the SZ angular power 
spectrum is proportional to the product of the MF evaluated at specific points 
in the mass-redshift space and the temperature (through the y-parameter de- 
pendence). Suppose that the normalization parameter erg is extracted by fitting 
a theoretical MF to an observed temperature function using a TMR Ti(M), 
and the gas temperature used in the evaluation of the cluster Comptoniza- 
tion parameter is derived from a TMR T2(M), such that Ti(M) > T2(M). As 
seen before, the Ti(M) scaling predicts lower cluster abundances than would 
be predicted employing the T2(M) scaling, and the T2(M) scaling attributes 
lower temperatures to a given mass. This could lead to an appreciable under- 
estimation of the SZ CfS. In the opposite case, when T2(M) is used to extract 
as, and Ti(Af) to evaluate the temperature, Ce will likely be overestimated. 
This inconsistency underlines the need for a unique choice of the TMR. 

2.4 Angular Power Spectrum 

The angular power spectrum of the SZ effect {e.g. , Molnar & Birkinshaw 2000) 
in the cluster population is characterized by 

j^^J-J ^(^' ^) ^) dM dz, (11) 

z M 

where N{M, z) is the MF, and r is the comoving radial distance. The function 
C(M, z) is the angular Fourier transform of the spatial profile of AT, which in 
the nonrelativistic limit is given in terms of the profile of the Comptonization 
parameter (cq. 5). In our calculations we consider clusters in the mass range of 
10^^ — 10^^ Mq h"^, but discuss how the results change when a narrower range 
is taken. The redshift integration interval - determined by the earliest epoch 
when clusters can be assumed to have formed and virialized - is uncertain and 
clearly depends on the cosmological model. 



3 Results and Discussion 

The 'standard' model in our calculations of the SZ angular power spectrum 
is the currently favored fiat, vacuum-dominated universe. The model is char- 
acterized by the following set of cosmological parameters: Jl^ = 0.3, a = 
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0.7, ^Ib = 0.05, n = l,h = 0.7. Clusters are assumed to be distributed accord- 
ing to the PS MF, have isothermal gas over a spherical region with a radius of 
lOrc, with the gas mass fraction that scales with the cluster mass and evolves 
in time as /(M, z) oc M^-'^ t{zy-'^^. We have also used two other recently pro- 
posed MFs, and generalized the gas distribution to a polytrope with a range 
of values for the index. In the calculations we have also considered the case of 
no gas evolution. The various cases of gas parameters for a given cluster are 
treated self-consistcntly by keeping the total mass and the gas mass fraction 
fixed as the polytropic index is varied. We carried out the calculations both 
without and with the SZ spectral dependence (the latter only for our standard 
model) evaluated at five observing frequencies - taken to correspond to the 
four MITO channels (De Petris et al. 2002) and the 545 GHz channel of the 
HFI on the Planck satellite. 

As mentioned in section 2.3, a mass-temperature scaling is needed in order to 
fit a MF to the observed cluster X-ray temperature function, so as to obtain 
the normalization as- We have adopted the mass-temperature relation and the 
deduced values of erg obtained by Wu (2000) from fits of the LS and ST MF 
to the observed temperature function derived by Viana & Liddle (1999). The 
latter authors relaxed the assumption of spherical gravitational collapse, and 
derived analytic expressions within the framework of non-spherical collapse. 
(We note that in these models a halo may begin to form when an overden- 
sity collapses along one of the principal axes of the deformation tensor. The 
earlier formation of clusters would delay the active formation epoch in com- 
parison with the PS case [Wu 2000], i.e., a cluster with present virial mass 
M has a higher probability of forming at a lower redshift than in the PS 
models.) For our choice of the standard model parameters, the values of as 
are 0.974,0.839,0.913 for the PS, LS and ST MFs, respectively. This param- 
eter was recently deduced also from the SDSS and WMAP CMB anisotropy 
databases, yielding somewhat lower values than we have adopted here (Mel- 
chiorri ct al. 2002, Bennett et al. 2003); however, the range of uncertainties 
in the various deduced values of as is wide and certainly includes the above 
three values that wc use here. Results for the power spectrum for the standard 
model with full, partial, and no gas evolution (to be specified below) are plot- 
ted in figure 1. Note that here and in the following figures the power spectrum 
is expressed in its dimensionless form, namely it is AT/Tq (rather than AT") 
that is expanded in Fourier space. In Figure 2, power spectra are shown for two 
polytropic models with indices 7 = 1.1, 1.3, and for the standard isothermal 
model. The predictions from the PS, LS and ST mass functions are shown in 
figure 3. The primary and SZ induced anisotropy power spectra are shown in 
figure 4. 

The SZ power spectrum in the standard model and in the models with the LS 
and ST MFs are shown in figure 3. Note that all three curves lie close to each 
other, particularly the PS and ST MFs. The LS MF yields lower values for 
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the power spectrum; as can be seen in the right panel of figure 3, this amounts 
to a reduction of the power spctrum by ~ 50%, somewhat higher at low mul- 
tipoles, and lower by ~ 10% at low multipoles; this curve virtually coincides 
with the PS prediction at higher multipoles. Apparently, when these MPs are 
fit to the same observational data by using the same mass-temperature rela- 
tion, differences among the power spectra will be minor. This is in particular 
important due to the fact that the MP is very sensitive to the input param- 
eters, in particular to those appearing in its exponential part. A substantial 
amount of effort has been made over recent years to improve the PS model 
by resolving some of its problematic features, and to obtaining more realistic 
MPs. It appears that at least for the task of calculating SZ power spectra, 
using these functions instead of the PS MP does not affect much the results, 
particularly so when compared with other sources of uncertainties. 

In all the models considered here the SZ power spectrum exhibits the same 
general behavior: a monotonic increase from low levels on large scales (low £), 
with a steep rise to a peak centered around I ~ 1500 — 3000 (depending on the 
model), followed by a steep decline. Differences between the respective curves 
are manifested in the orientation of the curve relative to the multipole axis 
and the magnitude - or height - reached at the peak. This characteristic curve 
stems from the angular Pourier transform of the Comptonization parameter. 
Por small values of I the Pourier transformed profile is essentially constant in 
magnitude; the parabolic shape is due to multiplication by £(£ + 1). On smaller 
scales the Pourier transformed profile is exponentially damped, becoming sen- 
sitive to the internal structure of the cluster. This gives rise to the steep fall of 
the power spectrum; the peak of the curve roughly represents the transition 
between large and small scales, or between unresolved and resolved clusters. 
Note that the power spectrum is calculated only for the thermal component 
of the SZ effect; the related anisotropy induced by the kinematic component 
has been treated in various papers {e.g., Cooray 2001). The magnitude of 
the integrated contribution of this component is second order in the cluster 
velocity (in the CMB frame), and thus very much smaller than the respective 
contribution of the main, thermal component. 

Gas evolution is a major factor determining the shape of the power spectrum, 
as is apparent from figure 1; with our adopted form of gas evolution the curve 
peaks at £ ~ 1400, reaching a magnitude of ~ 2.2-10"^^. Without gas evolution 
the peak shifts considerably to £ ~ 8000, where it reaches the much higher 
value of ~ 2.4 • 10~^°. Comparison of the two curves reveals that the higher 
the multipole, the greater is the contrast between the powers of the non- 
evolving and evolving gas. This stems from the higher gas densities in distant 
clusters predicted in the non-evolutionary gas scenario, which gives rise to 
an increase of power on smaller scales, reflecting the small angular sizes of 
distant clusters. It is also interesting to calculate the power spectrum in two 
'partial' gas-evolving models, the first chosen to correspond to ^ = and 
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r) = 0.2, and the second to ,^ = 1.45 and 77 = - namely removing from 
the gas mass fraction the redshift and mass dependences, respectively - in 
order to gauge the effect of the range of values of these parameters on the 
calculated power. For example, removing the redshift dependence results in a 
considerable dependence of the results on the mass lower limit, since - as can 
be seen in figure 1 - the non-evolving gas model still dominates the partially 
evolving gas model by a factor of ~ 17 at the highest multipoles. This maximal 
factor corresponds to the lowest cluster mass considered in the calculation, 
M — IO^^Mq h^^; this factor would have been smaller had we taken a higher 
mass cutoff. 

The f] = curve provides insight into the redshift range where most contribut- 
ing to the anisotropy. Since the power spectrum in the model with non-evolving 
gas is only a factor of ~ 2 larger (on average) than that of evolving gas, we can 
evaluate the maximal redshift at which the contribution is still appreciable: 
In a flat universe we have 



In 

t 




</2(l+z)3/2 
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(12) 



Now, the ratio between the power spectra excluding gas evolution and in- 
cluding only its temporal dependence is, {t/to)~^, and with the best fit value 
= 1.45, we determine the maximal factor this term contributes to the power 
spectrum in the non-evolving scenario to be (t/to)-^-^^ ~ Employing equa- 
tion (12) to solve for the redshift yields z ~ 0.24. The relatively low value of 
this redshift may seem surprising as we may expect high-redshift objects to 
contribute considerably to the power spectrum when there is no gas evolution. 
Instead, we only see a slight enhancement of its intensity, which peaks around 
a rather low value. One may then conclude that the ever decreasing cluster 
population at higher redshifts is responsible for this effect, and that the scaling 
of IC gas with mass dominates over its temporal change. 

It is unrealistic to expect that thermal conductivity in the gas is sufficiently 
high for it to be isothermal over the full extent of the gas distribution. A more 
general distribution is that of polytropic gas a polytropic equation of state, 
P oc p^, for which the Comptonization parameter y scales as a los integral of 
rf. Clearly, for isothermal gas and P — 2/3, 6 — 1 (eq. 2). We consider two 
values of the polytropic index: 7 = 1.1, 1.3, in addition to the isothermal case. 
Using the hydrostatic equation, the central temperature is 
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the functions ilj'{p,j3) and ou{p,f3) are defined in Colafranccsco ct al. (1994). 
Wfiile the central temperature is not sensitive to the exact value of 7, the 
Comptonization parameter has a steeper distribution in the polytropic case, 
resulting in a substantially different power spectrum than that of isothermal 
gas. Figure 2 describes the SZ angular power spectrum in the isothermal and 
the two other polytropic cases. While there is a slight increase of the power 
for 7 = 1.1 (as compared to the isothermal case), the peak value for 7 = 1.3 
is a factor ~ 2 lower with respect to the isothermal gas. The peaks shift 
to higher values of i in both ploytropic cases. These changes result from a 
steeper Comptonization parameter due to a more compact cluster, leading 
to smaller cluster size and thus an increase of power on smaller scales. The 
peaks of the two curves for polytropic models with 7 > 1 occur in the range 
£ « 2000—3000; the steeply falling part of the profiles shifts further to the right 
along the multipole axis. The difference between the magnitudes in the two 
polytropic cases is likely due to an interplay between two competing effects: 
the steeper Comptonization parameter assumes higher values in the central 
cluster regions, but it assumes these values along shorter distances. Therefore, 
with polytropic indices slightly above unity, there is an increase of magnitude, 
where the steeper profile dominates over the somewhat shorter distances along 
lines of sight with constant values of y, while for higher values of 7 the impact 
of the reduced pathlength dominates. 

In figure 5 we plot the primary anisotropy, calculated using the CMBFAST 
code, and show the results from WMAP, CBI, and ACBAR measurements, 
as well as SZ power spectra in the isothermal and polytropic cases considered 
above. We have included the spectral dependence which we calculated at a 
frequency of 31 GHz, as appropriate for the CBI experiment. ACBAR mea- 
surements (Kuo et al. 2002) were at a frequency of 150 GHz; the data points 
are shown as reported (in the latter paper). Were it known that the signal 
measured at high multipoles is largely due to the SZ effect, we could have eas- 
ily scaled the ACBAR data to the lower CBI frequency by the theoretically 
predicted spectral factor of ~ 4.2 (assuming a mean gas temperature of 7 
keV). This would have resulted in an enhancement of the ACBAR measured 
power over the data collected at the actual observing frequency. 

The 'excess' power detected by the CBI experiment at high multipoles (Mason 
et al. 2002) may possibly be attributed to SZ signal. It has been suggested 
that this measured 'excess' is likely to require a modification of as, in order to 
increase the cluster population, as predicted by the PS mass function. How- 
ever, such a modification is not that simple to make; as is usually estimated 
by fitting a theoretical mass function to the observed temperature function of 
X-ray clusters; therefore, its value is set by the fit. Moreover, different analy- 
ses yielded various values which span the rather broad range, 0.7 ^ o"8 < 1. 
An arbitrary modification of this parameter would have consequences on the 
mass-temperature relation, which might be in total disccordance with the 



13 



original mass-temperature relation employed in order to extract as in the first 
place. Our results seem to be consistent with the CBI measurements, thus 
circumventing the need for a modification in this parameter. However, even 
if this was not the case, discrepancies between the two sets of results may be 
explained by other means than the normalization parameter as- In particular, 
the not fully known properties of IC gas provide a certain degree of freedom 
in adjusting the theoretically predicted level of the anisotropy. For exmaple, 
as we have shown above, the anisotropy depends apprecaibly on the evolu- 
tion of IC gas; moderate (or no) evolution results in a considerable increase 
of the power spectrum at high multipoles. This is particularly relevant given 
the observational indications that there is weak or no evolution of the X-ray 
luminosity function of clusters. 

In figure 4 we compare the SZ power spectra obtained in the ACDM model 
with the corresponding power spectra of the primary anisotropy. The latter 
was computed using the CMBFAST code (Zaldarriaga & Seljak 1996). As can 
be seen from the figure, the SZ power spectrum begins to dominate over the 
primary anisotropy at multipoles of £ ~ 2000, but even at lower multipoles 
{£ > 1500) the magnitude of the SZ power already constitutes a fraction 
(> 10%) of the primary power. Furthermore, these results were obtained with 
constant spectral function, s{x) = 1. Its value in the Rayleigh- Jeans region is 
—2, but it can be much larger; for example at the highest frequency of the HFI 
instrument of Planck (~ 857 GHz), s{x) ~ 11, so - from merely a theoretical 
point of view the - there is a clear advantage in measuring the power at high 
frequencies. 

The power spectrum was calculated at five frequencies, four of which corre- 
spond to the channels, x = 2.5, 3.9, 4.8, 6.2, and a high frequency - the highest 
HFI channel with reasonably adequate sensitivity to SZ measurements aboard 
the PLANCK satellite, u — 545 GHz, x — 9.6. To assess the degree of impor- 
tance of using the relativistic SZ expression we plot the SZ angular power 
spectrum obtained using the non-relativistic expression, the exact relativis- 
tic expression, and the approximate analytic expression of Itoh et al. (1998): 
figure 6 compares the SZ angular power spectrum calculated with the non- 
relativistic expression of the spectral distortion with that calculated using 
the exact relativistic expression. The results demonstrate (again) the need for 
using the exact relativistic expression. The approximate analytic expression 
is of limited use near the the crossover frequency, as shown in figure 7. For 
example, at x = 2.511 and x = 3.907 and high values of £, it yields values 
that are 25% — 200% larger than those calculated using the exact relativistic 
treatment. 

Finally, to determine the redshift distribution of the SZ power spectrum con- 
tribution of clusters, we divide the redshift integration range (0.01 < z < 6) 
into six intervals, and plot in figure 8 the fractional contribution of each in- 
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terval to the overall power spectrum at x = 2.5 (the corresponding plots at 
the other frequencies do not differ substantially). As can be seen in the figure, 
the contribution to the power spectrum from redshifts 2; > 1 is negligible at 
multipoles i < 1000. The contribution of clusters at 1 < z < 2 increases from 
~ 10% at ^ ~ 2000, to as much as 40% at the highest multipoles. Clusters at 
2 < z < 3 contribute only at the highest multipoles at a level < 7%, and the 
contribution of clusters at still higher redshifts is < 1%, which is the reason 
for its absence in the figure. These results agree with the analysis done based 
on the partial gas-evolution scenario discussed earlier in this section. These 
results also imply that even though high gas temperatures are predicted at 
2; > 1 by the virial mass-temperature relation, which is based on a naive red- 
shift scaling, the impact of this unrealistic relation on our results is minimal 
since the contribution of high temperature clusters is negligible. 



4 Conclusion 

One of the main goals of our work has been to understand the reasons for the 
substantial variance among previous predictions for the power spectrum of 
the SZ-induced anisotropy. Among the various works on the power spectrum 
which employ analytic forms of the mass function to calculate the power spec- 
trum in the context of a given cosmological and DM model, the main reason 
for quantitative differences was shown to be a choice of the cluster mass- 
temperature relation that is inconsistent with value of erg which is adopted 
in the calculation. We have taken explicit account of the evolution of the gas 
fraction in the context of a simple model, and explored how the power spec- 
trum depends on the index of a polytropic equation of state. Our results are 
therefore self-consistent, more general in the descriptionof the properties of the 
gas than previous treatments, and include detailed predictions of the power 
spectrum at five frequencies based on a relativistic calculation of the spectral 
distribution of the SZ effect. 
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Fig. 1. The SZ angular power spectrum in the ACDM model, with gas evolution 
(upper left panel) and compared to a non-evolving gas scenario (upper right panel) . 
The lower left panel depicts the ratio between the two power spectra; the lower right 
panel describes the ratio between the power spectrum obtained in the non-evolving 
IC gas scenario and two 'partial' IC gas evolving models. 
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Fig. 2. SZ angular power spectrum for isothermal {5 = 1.0) gas and two polytropic 
temperature profiles with S = 1.1 and 1.3. 
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Fig. 3. SZ power spectrum calculated using the PS, LS and ST mass functions 
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Fig. 4. SZ and primary anisotropy power spectra in the ACDM models 
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Fig. 5. Recent results from the CBI, ACBAR, and WMAP experiments, together 
with the angular power spectrum of the primary anisotropy as calculated by CMB- 
FAST, and the SZ power spectrum in isothermal and polytropic models. WMAP, 
ACBAR, and CBI data are marked by dots, crosses, and empty circles, respectively. 
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Fig. 6. Power spectra calculated using the non-relativistic and exact relativistic SZ 
intensity change are shown in the left panel. The spectral dependence was calcu- 
lated at five frequencies, x = 2.5, 3.9, 4.8, 6.2, 9.6; each pair of adjacent solid and 
dash-dotted lines correspond to the relativistic and non-relativistic calculation, re- 
spectively. Line thickness increases with the value of x. Ratios of relativistic to 
non-relativistic power at each frequency are shown in the right panel. 
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Fig. 7. Comparison of the results from the exact relativistic and approximate an- 
alytic expressions. Solid lines show the power spectrum calculated using the exact 
relativistic calculation of Rephaeli (1995b) for the spectral function, whereas the 
dash-dotted lines show the corresponding results when the analytic expression of 
Itoh et al (1998) is used. Ratios of the results from the exact relativistic calculation 
to those from the analytic approximation are shown in the right panel; line thickness 
corresponds to frequency as in the previous figure. 
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x=2.511 




Fig. 8. Relative contribution to the overall power spectrum, at x = 2.5, by clusters 
at three redshift intervals. 
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